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Abstract 

In this paper, we study the realizabihty problem for retarded functional differential 
equations near an equilibrium point undergoing a nonlinear mode interaction between 
a saddle-node bifurcation and a non-resonant multiple Hopf bifurcation. In contrast 
to the case of transcritical/multiple Hopf interaction which was studied in an earlier 
paper the analysis here is complicated by the presence of a nilpotency which in- 
troduces a non-compact component in the symmetry group of the normal form. We 
present a framework to analyse the realizability problem in this non-semisimple case 
which exploits to a large extent our previous results for the realizability problem in 
the semisimple case. Apart from providing a solution to the problem of interest in 
this paper, it is believed that the approach used here could potentially be adapted to 
the study of the realizability problem for toroidal normal forms in the general case of 
repeated eigenvalues with Jordan blocks. 
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1 Introduction 



Retarded functional differential equations (RFDEs) are frequently used as models for various 
phenomena Pfl |21 CHI UH UHl UH 1201 • While the phase space for the resulting dynamical 
system is infinite dimensional, the existence of finite-dimensional invariant center manifolds 
near bifurcation points imply that much of the machinery developed for the analysis of finite 
codimension bifurcations in ordinary differential equations (e.g. normal forms, unfolding 
theory) are portable to RFDEs. Indeed, there is ample evidence in the literature, e.g. 
nBlEl! that these tools and techniques of local analysis can give much valuable information 
about the dynamics of RFDEs. 

In this context, one of the fundamental questions concerns the characterization of the 
range of dynamics accessible near a bifurcation point for a given RFDE model. This is 
not a trivial question, since even for scalar RFDEs, it is the possible to have bifurcations 
of equilibria with large dimensional center manifolds. In this case, if there are not enough 
independent delay terms in the nonlinear part of the RFDE, there may be severe restrictions 
on the possible dynamics which can be realized in the center manifold equations. The study 
of these types of questions is known generally as the realizability problem for RFDEs, and 
we refer the reader to the Introduction of |3] and to [HI IHl more details. 

In a previous paper we studied the realizability problem for scalar RFDEs in two 
important cases: the case of multiple non-resonant Hopf bifurcation, and the case of the 
interaction between a transcritical steady-state bifurcation and a multiple non-resonant Hopf 
bifurcation. In particular, we used the fact that these bifurcations admit normal forms with 
toroidal symmetries, and that these symmetries allow for a decoupling of the center manifold 
normal form equations into a radial part and an angular part. The uncoupled radial equations 
are a crucial part of the dynamics near the bifurcation, and it is thus very reasonable to 
consider the problem of realizability of the class of radial equations within a given RFDE. 
This problem was solved in general in ^ for the two important bifurcation scenarios described 
above. Moreover, it was shown in [3] that our solution to this realizability problem, that is, 
our estimate on the number of independent delays sufficient to achieve complete realizability, 
is optimal. The results in |3] considerably generalize previous results of Faria and Magalhaes 
jni and of Buono and Belair [3] . 

However, in Q the case of the interaction between a saddle-node steady-state bifurcation 
and a multiple non-resonant Hopf bifurcation was not studied, because there is a nilpotency 
associated to this bifurcation which is such that the approach used in |5 is not applicable. 
This nilpotency arises from the fact that the normal form for the saddle-node bifurcation 
contains an unfolding parameter which appears as an affine linear perturbation of the singular 
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vector field, i.e. 



P 



V 



0. 



A consequence of this nilpotency is that the normal form for the saddle-node/ multiple Hopf 
interaction admits a symmetry group which is a direct product of a compact torus group, 
and a non-compact group isomorphic to the additive group of real numbers. This latter non- 
compact portion does not occur in the cases studied in It thus becomes important to 
develop a framework which allows for the characterization of the effects of this non-compact 
symmetry on the normal form before being able to address the realizability problem for this 
bifurcation in RFDEs. 

In this paper, we adapt the Faria and Magalhaes normal form procedure [Ul [7j to the 
saddle-node/multiple Hopf interaction in scalar RFDEs, and use it to develop a framework 
suitable to studying the realizability problem for the radial part of the normal form in this 
nilpotent (non-semisimple) case. The framework is carefully constructed so that in the end, 
we may use our solution to the realizability problem in the semisimple case in jl] to the 
fullest possible extent. 

2 Preliminaries 

In this section we will briefiy recall some standard results and terminology in the bifurcation 
theory of RFDEs in order to establish the notation. For more details, see jH 13 113 • 

2.1 Infinite dimensional parameterized ODE 

Suppose r > is a given real number, n > 1 is a given integer and = C ([— r, 0] , M") is 
the Banach space of continuous functions from [— r, 0] into M" with supremum norm. We 
define G C„ as (^) = w (f + ^) , — r < 6^ < 0. Let us consider the following parameterized 
family of scalar (n = 1) nonlinear retarded functional differential equations 



from C\ into M and F is a smooth function from C\ x ]R'^+^ into M. In the "prequel" |3] to 
this paper, we assumed that F(0, 0) = 0, and DF{0, 0) = 0. While this hypothesis included 
the cases of multiple Hopf bifurcation and of a mode-interaction between a multiple Hopf 
and a transcritical type steady-state bifurcation, it excluded the case of a mode interaction 
in which the steady-state was of saddle-node type. Therefore, in this paper, we assume the 
following weaker hypothesis 



where L : Ci x 
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Hypothesis 2.1 F(0,0) = 0, DiF(0,0) = 0, and D^F(0,0) ^ 0. 

Performing a linear change of parameters and relabeling the parameters if necessary, we may 
then rewrite (|2.ip as 

z{t) = LoZt + u + F{zt,u,fx), (2.2) 

where we have set a = [u, fj.) G x W, Lq = L{0), and F{zt.,v,ii) = {L{u, /j.) — Lo)zt + 
F{zt, z/, fi) - DyF{Q, 0, 0). It follows from Hypothesis O that 

F(0,0,0)=0, DF(0,0,0) = 0. 

Clearly, the parameter v plays a distinguished role in ()2.2j) in comparison to the other s 
parameters 

Spectral hypothesis 

Suppose we set /i) = (0, 0) in (j2.2p . If we then write Lq as 

\Tl{d)(t){9), (2.3) 

-r 

where ?7 is a real- valued function of bounded variation in [— r, 0] and we let Aq be the 
infinitesimal generator of the semi-fiow associated with the linear RFDE z{t) = LqZi, then 
it is well-known that the spectrum cr(ylo) of Aq is equal to the point spectrum of Aq, and 
A G (tIAq) if and only if A satisfies the characteristic equation 



detA(A) = 0, A{X) = X- j dr]{e)e^^. (2.4) 

—r 

Denote by Aq the set of eigenvalues of a{Ao) with zero real part. 

Hypothesis 2.2 Throughout the rest of the paper, we assume the following hypotheses on 
Aq. Each element of Aq is a simple eigenvalue of Aq, and Aq has the following form: 

Aq = {0,±iuji,. . . ,±iujp}, 

where Ui, . . . , ujp, are independent over the rationals, i.e. if ri, . . . ,rp are rational numbers 
such that Yl^j=i ^j^j ~ 0' then ri = ■ ■ ■ = rp = 0. We further assume that the rest of the 
spectrum of Aq is hounded away from the imaginary axis. 
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Phase space decomposition 

In order to properly analyse the role of the parameters in equation ()2.2j) we need to 
augment this equation by considering the following system 

)= I' ^»^' + ''("'+-^ ((;;:)■'") I (2.5a) 

flit) = 0, (2.5b) 

where J^{{zt,iyt)'^ , fi) = F{zt,h'{0), fi{0)). The above system (j2.5|) can thus be viewed as an 
s-dimensional family (parameterized by /i) of 2- dimensional RFDEs. 

Taking into account Hypothesis 123 the linearized equation (i(t), i^(t)) = {L^zt + z^(0), 0) 
associated to ()2.5a|) has simple non-resonant characteristic values ± icui, . . . , ± iwp, and a 
characteristic value at of multiplicity 2. We then let P C C2 designate the 2p + 2- 
dimensional center subspace which is spanned by the columns of the following matrix 

... ij- ^2.6) 

We note that $ satisfies the linear differential equation — = where B is the {2p + 

dO 

2) X {2p + 2) matrix B = diag(0, iui, —iiOi, . . . , iup, —iup, 0) + N, and N is the nilpotent 
matrix whose entries are all zero except the entry at the intersection of the first row and last 
column, whose value is 1. 
We decompose C2 as 

C2 = P © Q. (2.7) 

Defining C2 = C([0, r], M^*), where M^* is the 2-dimensional space of row vectors, we intro- 
duce the adjoint bilinear form on x C2: 

(^, 0) = ^(0)0(0) - r r ^(^ - 9)dmmdC, (2.8) 



r JO 



where dr]{6) is of the form 
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drj^O) is as in (|2.3p and 5 (9) dO is such that 




We may then choose a basis {"iAi, • • • , '?/'2p+2} of the dual space P* such that if ^ = 
col('0i, . . . , V'2p+2) then <l>) = /2p+2, where throughout the paper we use the convention 
that for a given integer g > 1, is the q x q identity matrix. The following result will be 
useful later. 

Lemma 2.3 Let $ he as in \2. b]) and = {'4'k,i{C)) ^ (2p + 2) x 2 matrix whose rows 
form a basis for P* , and such that , $) = /2p+2- Then 

^2p+2,i(0) = 0, and i)kM ^ 0, 1 < A; < 2p + 1. 



Proof It follows from Hvpothesis 12.21 that 

f dr]{9) = 0, r dr]{9)e^'^^' = ±iu,, j = l,...,p 

J —r J —r 

and ^ ^ 

[ d7^ie)9^1, I dr]ie)9e^'-^'^l, j = l,...,p. 

J —r J —r 

Using ()2.8|1 . the equation $) = /2p+2 leads to 



^2p+2,l(0) 1 



1° dm = 0, 



^,,i(0)(^l-y dT^{e)ee''^^'j=l, A; = 2,4,...,2p, 

and ^ 

MO)(^l-J dr]{9)ee-'^^''^ =1, A; = 3,5,...,2p+1, 

so we get the desired conclusion. I 

Consider now the Banach space BC2 of functions from [— r, 0] into which are uniformly 
continuous on [— r, 0) with a jump discontinuity at 0. Elements of BC2 are written as 
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where G C2, A G and Xq is the 2x2 matrix- valued function 

( h = 

Xo{9) = { 

[ -r <e <0. 

Let TT : BC2 — > P denote the projection 

7r(0 + XoA) = <l>[(^,0) + ^(O)A], 
where ( , ) is the bihnear form ()2.8|) . We may then extend the sphtting ()2.7p to 

BC2 = P^keiTT, (2.9) 

with the property that Q C ker n. 

This structure now allows for a decomposition of the phase space which facilitates the 
implementation of the normal form procedure with parameters for RFDEs developed by Faria 
and Magalhaes in [7j. Specifically, it follows that (j2.5aj) is equivalent to the parameterized 
family 



5 ( ^ ) +^(0) ' 



X 



(2.10a) 



X 



(2.10b) 







where ^(0) is as in Lemma O x G R'^p+\ ly e R, e W, y e = Q n C\, {Cl is the 
subspace of Ci consisting of continuously differentiable functions), and Aq\ is the operator 
from Q} into kervr defined by 



A, 



+ ^0 



dm m - m 



2.2 Faria and Magalhaes normal form 

Consider the formal Taylor expansion of the nonlinearity terms in (|2.5p 



i>2 



where J^jiw) = Hj{w, . . . ,w), with Hj belonging to the space of continuous multihnear 
symmetric maps from C2+S x ■ ■ ■ x C2+S {j times) to M. If we denote x = {x, v)'^ and 
fj = ifjJj)^ where 



(2.11) 



\ 

/|(x,,,/.) = (/-.)Xo(-^^(^Y''''^ 
then ()2.10p can be written as 

i = 5x + 5^ ^//(x,y,//) (2.12a) 
i>2 ^' 

jy = AQ.y + Y.\f>^y^^) (2-12b) 

The spectral hypotheses we have specified in Hypothesis 12.21 are sufficient to conclude 
that the non-resonance condition of Faria and Magalhaes [7] holds. Consequently, using 
successively at each order j a near identity change of variables of the form 

(x,?/) = (x,y) + [/,(x,^) = (x,^) + (^/(x,^),[/2(x,/x)), (2.13) 

(where U^'"^ are homogeneous degree j polynomials in the indicated variables, with coefficients 
respectively in M^p+^ and Q^) system (j2.12|) can be put into formal normal form 

i = + ^ ^(7)(x,2/,//) (2.14a) 

j^y = AQiy + J2^g]{^,y,fi) (2.14b) 
i>2 ^' 

such that the center manifold is locally given by y = and the local flow of ()2.H1 on this 
center manifold is given by 

i = 5x + 5^ ^^7|(x,0,/i). (2.15) 
i>2 ^' 

The nonlinear terms gj in (j2.15p are in normal form in the classical sense with respect to 
the matrix B. 
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3 Non-semisimple Equivariant Normal Form 



The matrix B which appears in the previous section has a nilpotency associated with it which 
somewhat comphcates the computation of the normal form ()2.15|) . In particular, the analysis 
which was presented in j3] in the semisimple case does not carry over here. Nevertheless, 
in this section we will show how to generalize the normal form analysis presented in |lj to 
the present non-semisimple case. Apart from solving the problem which interests us in this 
paper, this approach may also shed some light on the general case where Hvpothesis 12.21 is 
relaxed to include repeated eigenvalues with Jordan blocks. 



3.1 X M normal forms 

Let ^E'(O) denote the {2p + 1) x 1 matrix obtained from the first 2p + 1 elements of the first 
column of ^(0) in Lemma (2.31 and let B be the {2p + 1) x {2p + 1) matrix 

B = diag(0, iuJi, —iuJi, . . . , iup, —iup). (3.1) 

It is easy to see from Lemma f2. 31 that fj in 1)2.111) is of the form 

/j(x,0,„)=(*(")^f--^) 

and we will thus consider normal forms for the following class of formal vector fields on 

/ x\ ( Bx + ubq \ ( fj{x,iy,fi) \ 

(3.2) 

where we will use the convention that is the unit row vector in R'' with all entries equal 
to zero except the entry in the row k + 1 whose value is 1, and where q will depend on the 
context. 

Using a mixture of complex and real coordinates, we identify 
]^2p+2+. ^ _ _ ^3,^^^^ jy^^^^ . . . ,//s) I Xo,v,iij e G C, j = 1, . . . ,s, A; = 1, . . . 

Define the following {2p + 2 + s) x {2p + 2 + s) matrix 

S = diag(5,0,), 

where 0^ is the s x s zero matrix, and let B^ denote the transpose of B. Let 
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where the closure is taken in the space of {2p + 2 + s) x [2p + 2 + s) matrices. Note that F 
is an abehan connected Lie group isomorphic to x R where is the p-torus: 

T^ = {diag(l,e^^^e-^^^...,e^^^e-^^^l,/,) | 9, e ^\ j = 1, . . . ,p} (3.3) 

and M is the one-parameter group parameterized as 

R = l2p+2+s + Q , OgM (3.4) 

where N"^ is the transpose of the nilpotent matrix 

N = B — diag(0, iui, —iui, . . . , iup, —iup, 0, O5). 

Definition 3.1 For a given integer i > 2 and a given normed space X, we denote by 
}j'^p+'^+^(^X) the linear space of homogeneous polynomials of degree t in the 2p-|-2 + s variables 
X = {xq,Xi,Xi, . . . ,Xp,x^), V and /i = (/ii, . . . ,/is), with coefficients in X. For X = M^P"'"^+'^, 
define i/f +^+"(M2p+2+^, r) C //f +^+"(M2p+2+s) ^/^g subspace of T-equivariant polyno- 

mials, i.e. 

/ G J/f+'+'(M2p+2+«,r) ^ 

/ e i/f+^+'(M2p+2+s) and lf{T^w) = f{w), = (x, z/, yu) G M^p+z+s^ V7 G T. 

Normal forms for (j3.2j) are computed using the homological operator 

jC-s ■■ H^/'^^-^\M.^P+^+') — > //f +'+'(M2p+2+^) 

/ ^ iCj)iw) = Df{w)Bw - Bf{w). 

To specify the normal form, we must find a complement in if^^^^^*(]R2P"'"2"'"*) to the range 
of Cj^. The following is a very well-known result in the theory of normal forms |^ IIUI ITT|. 

Proposition 3.2 

It is straightforward to compute the general element of if^^"^^^*(]R2p+2+*, F). 
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Lemma 3.3 A smooth vector field f : M2p+2+i5 — 
/ is of the form 



ai{xQ, XiXi, . . . , XpXp, jj) Xi 



2p+2+s T^-equivariant if and only if 



f{x,v,^) 



fl) Xp 

Qjpi^XQj XiXi^ • • • J ^p^pj ^5 /^) ^pi 

b X Q J X ]^ X ^ ^ • • • ^ X pX p J jy J ^ ^ 

Ci{xo,XiX]:, . . . ,XpX^, 



(3.5) 



y Cg (xq t XiXi, . . . , XpXp, ly, /i) , J 

where ai,...,ap are smooth and complex-valued, and Oq, 6, Ci, . . . are smooth and real- 
valued. Furthermore, a vector field of the form \3. ^) is T-equivariant (where F = x z/ 
and only if aQ,ai, . . . ,ap,ci, . . . ,Cs are v -independent, and 

o-o = go{xQ, XiXi, . . . , XpXp, fx), 

b = u go{xo,xix^, . . . ,XpX^, ^) gi{xo, xixi, . . . , XpX^, ^) 

for some smooth functions Qq and gi . 

Proof The vector field / is F-equivariant if and only if / is both T^-equivariant and 
M-equivariant. For the general form of the T^-equivariant vector field ()3.5|) . see |4.. We then 
further require that ()3.5p commute with all matrices of the form l2p+2+s + © N'^, O G M. 
The result follows after a straightforward computation. I 

The vector field / in ()3.2|) has the special form / = (/, 0,0), which we want our normal 
form changes of variables to preserve. Since we are only interested in the first 2p + 1 com- 
ponents of ()3.2|) . we would like to obtain a splitting of H'^^'^'^'^^ (R'^^^'^) akin to the splitting 
of i/^*'^^^'^(]R^P+^+*) in Proposition 13.21 For this purpose, we will need the following 
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Definition 3.4 



(a) We define H^/^^^' {R^p+\ Tp) to be the subset of Hf^'^^\R^'P+'^) consisting of mappings 
f : M2p+2+s — ^ ]R2p+i yjfiose components are of the form of the first 2p+l components 
of^. 

As the notation suggests, note that H'^p~^'^~^^{M.'^p~^^ ,Tp) consists precisely of elements 
0/ if|^^^^*(]R^^"'"^) which are equivariant under an action ofV on 

/ G //f+'+^(M2p+i) and /(70 a:, z/, /i) = jof{x,u,l^), V70 E Tq, V(x,z/,/i) G M^^, 

where Tq is the group of {2p + 1) x {2p + 1) matrices which is isomorphic to V , and 
is parameterized as 

^o = {diag(l,e^^^e-^^^...,e^^^e-^^'') | 6, e S\ j = 1, . . . ,p} (3.6) 

(b) We define i/f +^+'(R2p+i, T^) to be the subspace 0/ //f +^+"(R2p+i, T^) consisting of 
V -independent and -equivariant polynomials. 

(c) We define the following operator 

df 

f I — > {CB,u){f){x,u,fi) = D^f{x,u,ij,)Bx - Bf{x,u,ij,) + z/ — (x, z/, 

where B is as in and note that Cb,v is the usual homological operator associated 

to the X component of \3. Furthermore, we have Cj^{f, 0,0) = {Cb,u f, 0,0) . 

Proposition 3.5 

i7f+'+^(R2f+i) = //f+^+^(M2p+i, TP) ©rangers,,. 



Proof The proof is given in the appendix. 
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3.2 Equivariant projection 

We will now construct an appropriate linear projection associated with the splitting of 
^2p+2+s^^2p+i^ given in Proposition E^l 

Definition 3.6 Let / d'j denote the normalized Haar integral on Fq = V (see 113.6]) ). We 



To 

define the linear operator 



A : i/f +2+"(R2p+i) — > iff +2+"(R2p+i) 



/^(i/)(x,z/,/i)= / -ff{r'x,0,fi)dj. 



Proposition 3.7 A is a projection. Furthermore, 

range A = iff +'+'(M2p+i, jp) (3.7) 

and 



ker A = range Cb,u- (3i 



Proof The proof is given in the appendix. 

For any / G iif ^^^'(M^p+I), write 



f = Af + {I-A)f, 

and note that Af is T^-equivariant and //-independent and that (/ — A)/ G keryl. From 
Proposition Ha there exists h G iif ^^^'(M^p+I) such that £5,^/1 = (/ - A)f. 

3.3 Phase decoupling 

Elements of the space 

iif+'+^(M2p+i,TP) 

are i/-independent and equivariant with respect to 
the torus group Fq = defined in ()3.(i|l . As seen in [3], this toroidal equivariance can be 
used to achieve a decoupling of the normal form. Specifically, we have 

Proposition 3.8 Consider the following differential equation on M.'^f^^: 

X = B X + UBq + f{x,fi), 
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where f is smooth, v -independent, satisfies /(0,0) = 0, 15/(0,0) 
equivariant, i.e. f has the form 



0, and is Tq = JP- 



^ 1 ^ ^0 J 1 1 5 * * * 5 X pX p J ^ ^ X ]^ 



(3.9) 



Q'p('^O) 2^1X1, . . . , XpXp, fl) Xp 
y Q.p(xo5 XiXi, . . . , XpXp, /i) Xp, y 

where ai,...,ap are smooth and complex-valued, and is smooth and real-valued. Then 
under the under the change of variables xq = po, xj = pjc^^^ , j = 1, . . . ,p, this differential 
equation transforms into 



Po = z/ + ao(po,P?,---,Pp,At) 
Pj = Re{aj{po, pI,..., pf,, p)) pj, j = l,...,p 



and 



(3.10) 



(3.11) 



Oj = iujj + lm{aj{po,pl, . . . ,Pp,/i)), i = 1, . . . ,p. 

Proof This is a simple computation. I 

As in i4 , we will call the subsystem ()3.10|1 the uncoupled radial part of the normal form 
(I3.9|l . Recall that this uncoupled radial part has some residual reflectional symmetry. Denote 
by Z2,p the group whose action on R^^+^ is given by 



{po, Pi, • • • , Pp) (Po, Aipi, . . . , Appp), 
where \j E {1,-1}, j = l,...,p. 



(3.12) 



Definition 3.9 For a given integer £ > 2 and a given normed space X, we denote by 
^p+i+s^j^^ t/ie linear space of homogeneous polynomials of degree t in the p + 1 + s variables 
p = (po, pi, . . . , Pp) and p = {pi, . . . , ps) with coefficients in X. Denote by H^'^^^^iW''^'' , ^2,p) C 
^p+i+s^j^p+i^ i/ie subspace of H'^^'^'^^iW'^^) consisting of Z2^p- equivariant polynomials. 
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It is easy to show (see ^T]) that the most general element of H^^^^'^ {W^^ , ^2,p) has the form 
of the right-hand side of (|3.1(jp . We then define the following surjective linear mapping 

n : iff +'+'(M2p+i, TP) — , H^+^+'{RP+\ Zs.p) (3.13) 

which is defined by sending the general element ()3.9p of H'^^^^^^ [M.'^p^^ ,V) to the following 
element of H^+'^^'{W+\ Za.p): 



/ ao(po,P?, • • • ,Pp,/i) \ 
Re(ai(po,P?,---,Pp,/^))pi 

\ Re(ap(po,P?,---,Pp,^))Pp / 



(3.14) 



3.4 Parameter splitting 

As seen in it is useful when considering unfoldings to be able to refine Propositions 13.51 
and 13.71 in order to make explicit the roles of [x, u) as primary variables and n as unfolding 
parameters in ()2.14ap . For this purpose, we define the following spaces. 

Definition 3.10 Let d > 1 be a given integer (for our purposes, d will he equal to either 
2p+l or to p+1), £ > 2 be an integer, X be a normed linear space, and G be a linear group 
acting on X . Let Hf^^^'"^ (X) be the linear space of homogeneous polynomials of degree i in 
the variables {^i, ■ ■ ■ , ^d) , ^ o,nd fii, . . . , fXg, and Hf~^'^{X, G) C Hf~^^~^'^(X) be the subspace of 
1/ -independent and G-equivariant polynomials. 

(a) We define Hf+\X) C Hf+^+'{X) and Hf{X,G) C Hf+\X,G) to be the subspaces of 
^-independent polynomials. 

(b) We define Pf+^+'{X) C Hf+^+'{X) and Pf+'iX, G) C Hf+'{X, G) to be the subspaces 
of polynomials which vanish at fi = 0. 

It follows from these definitions that 

i/f +'+'(R2p+i, JP) = H^p^\R^p+\ TP) © +'+'(M2p+i, TP) (3.15) 

HP+'+'(RP+\Z2,p) = Hl^\W+\Z2,^®P!^'^\W+\Z^^,). 

Furthermore, the various operators Cb,u, A and 11 previously defined preserve these decom- 
positions. We then get the following refinement of Propositions 13.51 and 13.71 
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Proposition 3.11 

n(//f +'(R2p+\ TP)) = //J+'(RP+i, Z2,p), n(P^'P+'+'(R2p+i, TP)) = P/+'+'(RP+i, Z2,p) 

Consequently, if we define A\i and A\2 to he respectively the restrictions of A on Hf'^'^ {E?^^^) 
and on Pf'^^'^^'^ (M^'p^^) , then A\i and A\2 are projections. Similarly define Cb,u\i and Cb,u\2 
to be respectively the restrictions of Cb,u on //^^^^(R^P"*"-^) and on P^^P'^^^''(R^P"'"^), then 

rangeili = H^p^\R^p+\ Jp), rangei^ = Pf +'+'(R2p+i, Jp)^ 

kerv4|i = range Cb,i^\i, kery4|2 = range iZ^.i^b, 

^2p+2^^2p+i) _ (M2p+i^ TP) ©range 

P/P+2+^ (R2P+1) = pfP+^+-\m^P+\TP)(B range Cb,A2- 

We now combine the results of this section with the Faria and Magalhaes normal form 
procedure outlined in Section 2 in order to obtain the following version of Theorem 5.8 of 
[H] and Theorem 2.16 of |7j. 

Theorem 3.12 Consider the system \2.1'^) 

d . '"4- .2, , (3.16) 

—,y = Aqiy + fA^,y,ii), 



dV 

i>2 



where x = ( ^ ) . Write 



f}Mf^) = ( ) = ( ^^("^ 7^-("'^) ) , (3.17) 

where hj G i/J^^^(R^P^^) and qj G P^^^^^^'*(R^p+^). Then there is a formal near-identity 
change of variables 

(x,y) -> (x,y) + iU\^),U\^)) + {W\^,fi),W\^,fi)) 



16 



(where VF^(x, 0) = and Vr^(x, 0) = 0) which transforms Vj.lt^) into system \2.14\} (upon 
dropping the hats ), and the flow on the invariant local center manifold y = is given by 

X = Sx + i^eo + ^((i|i(/ij + y,))(2;) + (^|2(gi + ^i))(2:,/x)), 

J>2 

z> = (3.18) 
/i = 0, 

where Y2 = 0, Z2 = 0, and for j > 3, Yj = Yj{x,iy) and Zj = Zj{x,u,fi) are the extra 
contributions to the terms of order j coming from the transformation of the lower order 
(< j) terms, and Zj vanishes at n = 0. 



4 Main results 

In the semisimple case of jl], the main reahzabihty results are a consequence of the surjec- 
tivity of a certain hnear operator between suitable spaces of polynomials, which is proven 
in Proposition 4.3 of |3]. In this section, we define the corresponding linear operator in the 
present non-semisimple case, and prove that its surjectivity follows from the surjectivity in 
the semisimple case. Once this has been achieved, we will state our main realizability results 
for the class of uncoupled radial equations ()3.14|) which are obtained from writing the normal 
form center manifold equations (j3.18|) in polar coordinates. 



4.1 Linear analysis 

Definition 4.1 For a given integer i > 2, let H'^^'^'^^^ (M.) denote the linear space of homo- 
geneous degree £ polynomials in the 2p + 2 variables 



V = (vo, vi, . . . , Vp 



Wo J ' \ Wi J ' " ' ' \ Wp 



and s parameters jj, = (/ii, . . . , /i^) with real coefficients. Denote by H^^^^^{M.) the linear 
space of homogeneous degree i polynomials in the variables v = {vo,vi, . . . ,Vp) and fj, = 
(/ii, . . . , /is). Note that we may decompose these spaces as in i3.1^} as follows 

i/f +'+^(R) = /Jf +'(M) © +'+^(M) 

and 

(R) = HP+\R) © P/+'+'(R) 
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where H^p+\R) and Hf+\R) are the ^-independent polynomials , and P^^'^'^^'^ (M.) and Pf'^^^'^ (R.) 
are the polynomials which vanish at = 0. Finally, define the surjective linear mapping 



as 



{'JZ{h)){v,^) = h 











(4.1) 



Let $ be as in fj2.6|) . We will define $ to be the 1 x {2p + 1) matrix obtained from the 
first 2p + 1 elements of the first row of $, i.e. 



From Lemma f2. 31 and the fact that ^1^(0) denotes the {2p + 1) x 1 matrix obtained from the 
first 2p + 1 elements of the first column of ^(0) in Lemma f2. 31 it follows that 



Vl>(0) = Col(Wo, Ui,Ui, . . . ,Up, Up) 



(4.2) 



where 7^ is real and Uj 7^ are complex, j = 1 
Let r = (ti, . . . , Tp+i) G M^"*"^, and define 




( 1 



























1 


1 




g-iaJir2 






^2 

















1 


1 








g-ia;pTp+i 


















1 



(4.3) 



]_ gi<^lT2 g-«<^lT2 



gioJpTl g-iWpTl \ 

^iuJpT2 g — «'^pT2 



v 



(4.4) 



gltJlTp+l g — liJlTp+1 _ _ _ gJOJpTp+l g — JOJpTp+l 

As in we define the ^-mappings associated to Er and to E-j- respectively as follows: 

(^^(/i))(x,/i)^vl/(0)/i(i,x,/i) 
18 



(4.5) 



and 

(4.6) 

Now, let n : //f +^+'(M2p+i^ -jp) — > H^^'^-^'{W+\Z2,p) be the mapping defined in (jSISD- 
dSH, let A : i/J^+^+'(M2p+i) — , //f+^+"(M2p+i,TP) be the projection operator defined in 
Definition 13.61 and define 

by 

A9{x,fi)= 75f(7"^x,/x)ci7, 
as in (4j. Our main result in this section is the following: 

Proposition 4.2 For an open and dense set U C R^^"*^, the following linear mapping is 
surjective for all t ElA: 



Proof It is a simple computation to show that 

n o i o = (n o A o o 7^, 

where TZ is the surjective linear operator defined in ()4.1|) . In it was shown that there 
exists an open and dense set U C W'^^ such that for all r eU, the mapping 

is surjective. Thus, for all r G W, 11 o A o is surjective. ■ 
4.2 Main results on realizability 

With Proposition 14.21 in hand, we obtain realizability results analogous to the semisimple 
case. Since the proofs are almost identical, we merely give the statements of these results 
here. First, we will define the following linear spaces of (non-homogeneous) polynomials 
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Definition 4.3 



Theorem 4.4 Consider the RFDE \2.^) . and let Aq denote the set of solutions of \2.4\) with 
zero real part. Suppose that Hypothesis \2.2\ is satisfied. Let i > 2 be a given integer. For 
each h G 7{^'^^ {W^^^ , Z2^p) and each q G V^^^^'^{W^^,Z2^p) there are p + 1 distinct points 
Ti, . . . , Tp+i G [-r, 0], an rj G and a ^ E Pf+^+'(M) such that if 

F{Zt, /i) = r]{z{t + Ti),...,z{t + Td)) + ^{z{t + Ti),. . . ,z{t + Td),n) 

in \2.^) . then in polar coordinates, the radial part of the center manifold equations VJ.lt^) in 
X 'R-equivariant normal form up to degree i reduces to p = v Bq + h{p) + q{p,p), where 
p = {po, pi, . . . , Pp) . In fact, T = (ti, . . . , Tp+i) can be chosen in an open and dense set of 
[—r,OY~^^, independently of the particular h and q to be realized (i.e. only rj and ^ must be 
changed in order to account for different jets to be realized). 

As in the semisimple case of j3] , we can show that the number of delays p+1 shown 
above to be sufficient to solve the realizability problem for the radial part is optimal i.e. 
surjectivity is violated beyond some finite order io if the number of delays is less than p+1 
(see Theorem 5.4 of i4j). 

For general RFDEs, we have the following result on realization of unfoldings: 

Theorem 4.5 Consider the general nonlinear RFDE 

z{t) = LoZt + u + N{zt) (4.7) 

where u E R, Lq : Ci ^ W is a bounded linear operator from Ci = C ([— r, 0] , M) into M, and 
N is a smooth function from Ci into M, with N{0) = 0, DN{0) = 0. Let Aq denote the set 
of solutions of ^2.4\ ) with zero real part and suppose that Hypothesis \2.2\ is satisfied. Then 
the local dynamics of j^. 7| j near the origin on an invariant center manifold can be described 
by a system of ordinary differential equations on M^^'^^. Moreover, this ODE system can 
be brought into x R-equivariant normal form to any desired order i, and the resulting 
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(truncated at order i) normal form can be uncoupled into an p + 1-dimensional system and 
a p-dimensional system 

p = ueo + h{p;N) (4.8) 
= k{p-N), (4.9) 

where for given N , h{- ; N) is some element of n^^'^\W+\ Z2,p) , and k{- ; N) : W — > W. 
Let h{p, p) be an s-parameter equivariant unfolding of h of degree at most i, i.e. h & 
/^p+i+s^]^p+i^ ^^^^-j and h{-,0) = h{- ; N). Then there exists an s-parameter unfolding of 
(|//. 7| ) of the form 

z{t) = Lo{zt) + + N{zt) + i{z{t + n), . . . , z{t + r^+i), ^) (4.10) 

(where r = (ri, . . . ,Tp+i) G M'''^^, and ^ G ^^^^^^(IR) vanishes at p = 0) which realizes the 
unfolded radial equations 

p = ueo + h{p,p) 
on an invariant center manifold for ^.1(J{ ). 

4.3 Generic saddle-node/Hopf interaction 

The following example is an immediate consequence of our results. 

Example 4.6 Consider the RFDE (j2.2|) in the case p = 0, 

z{t) = LoZt + u + F{zt), (4.11) 

such that the characteristic equation (j2.4|) has simple purely imaginary roots ±iu! 7^ 0, a 
simple root at 0, and no other roots on the imaginary axis. If 

F{zt) = A2o{z{t + n)f + Ai,z{t + n)z{t + r2) + Ao2{z{t + T2)f 

A3o{z{t + ri))3 + A2iz{{t + n)fz{t + T2) + Auz{t + n){z{t + ra))^ + Ao3{z{t + ra))^, 

(4.12) 

where ri,T2 G [— t, 0], then the uncoupled radial part of the center manifold equations to 
cubic order are the following Guckenheimer [T^ IT3] normal form 

po = /y + aipl + a2pi + aj,pl + a^popl 

(4-13) 

Pi = hipQpi + h2pl + hpipl, 
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where the coefficients ai, 2,3,4 and 61^2,3 are functions of ^420, An, Aq2, ^30, ^21, A12, Aq^, 
Ti and T2. From our results, it follows that generically, any values of ai, 2,3,4 and 61,2,3 can 
be achieved by appropriate choice of ^420, ^11, ^02, ^30, ^21, ^12, ^03- Also, the following 
versal unfolding of ()4.13p 

Pq = u + aipl + a2pl + a^pl + a^popl 

(4-14) 

Pi = PPi + hpoPi + hPi + hpipl, 
is generically realized by the following unfolding of ()4.11|) 

z{t) = Lozt + U + F{zt) + p {Aioz{t + Ti) + A01P2 z{t + T2)) 
for appropriate choice of Aiq and Aqi. 

5 Concluding remarks 

Our solution here to the realizability problem for the non-semisimple saddle-node /multiple 
Hopf interaction in scalar RFDEs complements our results in ^ for the semisimple cases of 
transcritical/multiple Hopf interaction and non-resonant multiple Hopf bifurcation. We note 
that our results here on realizability are generic, optimal in the number of delays required 
to guarantee realizability, and applicable to any finite order expansion and truncation of 
the normal form. Therefore, nonlinear degeneracies and their unfoldings for the saddle- 
node/multiple Hopf interaction are covered by our theory. 

As in [4, we have not considered general n > 1 dimensional systems of RFDEs. This 
is not based on any deep theoretical issues associated to the n > 1 case, but rather to 
complications arising out of notation and messy algebraic computations which would make 
the exposition extremely cumbersome. We have, however, every reason to believe that the 
realizability problem for each of the bifurcations studied in [HI and in the present paper could 
be solved for systems using the techniques and framework we have developed. 

More subtle is the issue of relaxing Hypothesis 12.21 to include repeated eigenvalues with 
Jordan blocks, and rational resonances in the purely imaginary eigenvalues. Apart from 
affecting the dimension of the torus group admitted by the normal form (and consequently the 
dimension of the uncoupled radial equations), Jordan blocks would introduce an additional 
non-compact component to the normal form symmetry. In this case, our present analysis 
in this paper may shed some valuable light on a suitable approach to tackle the general 
problem, i.e. ffist solve the semisimple case, and then try to exploit the semisimple solution 
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as much as possible in order to prove realizability in the associated non-semisimple problem. 
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A Proof of Proposition 13.5 



Let / be a given element of i/f +^+'(M2p+1), and consider / = (/,0,0) G //f +^+'(M2p+2+^). 
From Proposition 13.21 there exists h = {h^,h'^,h'^) G if^^^^^'*(M^P+^+*) and a unique g = 
(^^^/^^'') e fZ'f+^+"(R2p+2+«,r) such that 

(/, 0, 0) = r , hn + {g\ 9\ gn- (A.I) 

Note that 

/i^ h") = Cb,u /i^' - /i'' eo , D.^h'^Bx + v— , D^h^Bx + u 



\ OXq OXq 

Thus, from (jA.lll . we get 

D.^h'^^'^Bx + u— — + g"^" = (A.2) 
dxo 

and 

/ = /:B,./^"-/i"eo + ^^ (A.3) 
From Lemma [3. 3| g'^ is of the form 

g" = Z/ri(xo,XiXr, . . .,XpX^,fi) + r2(xo,XiXl, . . .,XpX^,fi), 

where ri is such that the component of g^ along eo, g^, is of the form 

9o ~ -^0 ''"i (^0) -^i-^i ) • • • 5 ^pXpi f^) ■ 
Therefore, the ly part of equation ()A.2|) reduces to 

dh" ^ . / dh" _dh"\ ^ ,^ ^, 

- + L [^^d^ - ^^d^ )+-n + r, = 0. (A.4) 

Before we proceed any further, we will need the following two lemmas. 
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Lemma A.l Let h be a real-valued smooth function of Xq, Xi,Xi, . . . , Xp, Xp, u and jj, such 
that the function 0(xo, xi, . . . , Xp, a^, z^, /i) defined by 



, , . / dh dh\ 

0(xo, xi, xi, . . . , Xp, Xp, z/, /i) = 2^ lUj I Xj- n= J 



(A.5) 



is invariant, i.e. 

0(xo, e^^^xi, e"*^ixr, . . . , e^^^Xp, e"*^^x^, i^, yu) = 0(xo, Xi, xT, . . . , Xp, x^, z/, /i) (A. 6) 

/or all 9i, . . . ,9p E^, and for all (x, z/, /x) G R^p+^+*. T/ien is also V invariant, and = 0. 

Proof A simple computation using ()A.5|) and ()A.6jl leads to 

Integrating this equation gives 

h{xo, e*'^i^xi, e"*'^i^xl, . . . , e"^f^Xp, e~*'^''^x^, u, /i) — /i(xo, xi, xT, . . . , Xp, x^, i^, /x) = 

^ 0(xo, xi, xT, . . . , Xp, x^, z^, /i) . 

For any given (x, ly, /i) G ]R^p+^+'^, the left hand side of the previous equation is bounded in 
^, and it follows that must equal and that 

The conclusion follows from the non-resonance condition on uoi, ... ,ujp specified in Hypothesis 
12 .21 using density and continuity. I 



Lemma A. 2 Equation II1A.4}) implies that must be of the form 

h" = [){xo,XiXi, . . . jXpX^jU, ^) = — ri(t, XixT, . . . , XpX^, /i) (it — 

Jo 

m(xiXi, . . . , XpX^, u, /i). 
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Proof Write h'^{x^v,^) = Yl^j=o cij{x,jjL)vK Equation ()A.4|) becomes 

where g = — (z/ri + is T*^ invariant as in Lemma lA.ll Applying Lemma lA.ll suc- 
cessively to the coefficient of , . . . , in (IA.7|) . we get that r2 = 0, 7—^ = — ri, 
dd ■ 

7-^ = 0, j = 1, . . . , and Oq, . . . , are T^-invariant. The conclusion of the lemma then 
follows immediately upon integration. I 



The component of ()A.3|) along Bq now has the form 
dhl >A , / dhf. _dh^\ 



fXo 



/ ri{t, xix^, . . . , XpX^, /i) dt + m(xix^, . . . , XpX^, + 
Jo 



(A.8) 



xo ri(xo,xiZi, . . .,XpXp,fi). 

Using Taylor's theorem, we write 

m{xiXi, . . . , XpX^, u, /i) = m(xiXi , . . . , XpX^, 0, /i) + m{xiXi, . . . , XpX^, i^, /x) , 

and note that ^ 

u m{xi'xi, . . . , XpX^, i^, /x) = z/— — xotn, and 



dxo 

_ d 

XjT^Xoxh ] = 0. 
ox] 



Consequently, if we denote /ig = /ig + xoxh, then ()A.8|) reduces to 

/o = '^^ + ^^uJj\x,-^-x--^] + 
^^0 ~l y (JXj ox] J 



I ri(t, Xix], . . . , XpO^, /i) dt + m(a;iaJi, . . . , x^a^, 0, + 



(A.9) 



?"i(xo, XiXi, . . . , XpXp, fl). 
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Together, ()A.9|) and the last 2p components of ()A.3|) imply that / can be written as a sum 
of an element in range £5,!^ and an element in Hf^^^^^{M.'^^^^ i.e. 

^2p+2+.^^2p+i) = iff+^+^(R2f+i,Tf) + range 
We now must show that the above sum is, in fact, a direct sum. Suppose that G 

CB,uh^ + ^7" = 0. (A.IO) 

Write the component of along eo as 

5fo(xo,a;ixI, . . .,XpX^,fi) = go{0,Xix^, . . .,XpX^,fi) + Xog^{xo, Xix^, . . .,XpX^,fi), (A.ll) 
and define 

'"l('^O) XiXi, . . . , XpXp, /i) 

' ' tg^{t, XiXi 



Jq \ ^ ^ dxo 

Note that ri is regular at = 0, since both the integral term above and its derivative with 
respect to xq vanish at = 0. A simple computation verifies that 



(A.12) 



■^0 9q{-^0j XiXi, . . . , XpXp, /^) = / f'lij': XiXi, . . . , XpXp, /i) dt-\- 

Jo 

xo ri(xo, xixI, . . . , XpX^, fi). 

Consequently, if we define g^ = 0, = 0, 

g" = iyri{xo,XiX^, . . .,XpX^,fi) 

and 

l>XO 

h" = - ri(t, XixT, . . . , XpX^, /i) dt - (0, XiX^, . . . , XpX^, /i), 
^0 

then from Lemma we see that {{xq ri,gf, . . . , g2p),g'', g^) £ i^f^*'^^^'^ (R^^+^+'* , F). It now 
follows from (jA.lOj) that 

(0, 0, 0) = /:^(/^^ /^^ /i^) + ((xo n, ^7?^, . . . , g^^^g^, g^). 
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From Proposition 13. 2( we get that ri = gl = ■ ■ ■ = g^p = Q and C j^{h^ , , h^^) = (0,0,0), 
from which it follows that 

i^B,u h"" + (0, XixT, XpX^, n) eo = 0. 
The component of the above equation along eg, is 

^ + '^■^^ [^^ d^~''^d^)^ ^ ' • • • ' 

Using Lemma lA.ll we conclude that g^lO /i) = 0, CB,vh^ = 0, and from 

()A.11|) . that g^ = 0. Therefore, we conclude that 



B Proof of Proposition 13.71 

The proof that ^4 is a projection is similar to the proof given for the projection operator A 
used in the semisimple case in p|. 

Now, let / G range A, then Af = /, i.e. 

f{x,u,fi)= f{'j-^x,0,fi)d-f. 

So, for any a G Fq, we have 

7 /(7~^x, 0, fi) d-f = f{x, V, fi) = f{x, 0, fi), 
and thus / G //f +^+'(M2p+i, tp). On the other hand, if / G H^p^'^\R'^p+\V), then 
{Af){x,v,y)= j 7/(7~^x, 0,^)(i7 = / f{x,0,n)dj = 

/(x,0,/i) = f{x,iy,fi), 
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so / G range A. This establishes (|3.7|) . We now estabhsh (|3.8|) . Since A is a projection, then 

^2p+2+s^^2p+i^ = range A © ker A. 

From Proposition 13.51 we conclude that dim ker A = dim range Cb,^- Thus, we need only 
show that range I, C ker A. Recall the following lemma which was proved in 

Lemma B.l Let g : Tq — > M.^p^^ be a continuous function, then 

I ^(7)^7= lim ^ ! g{e^')ds. 

Jto ^ Jo 

Now, let / G range £^^1,, then there exists g G i^f^^^^'''*(M^P"'"^) such that 

Og 

Therefore, using Lemma fB. 11 we get 

(i/)(x,i^,/i)= / 7/(7~'a;,0,^)(i7= hm 1 re^V(e-^'^,0,/x)ds 

1 

= lim - / e^^ iD^g{e-^'x,Q,ii)Be-^'x - Bg{e-^'x,Q,n)) ds 
= \im ^ f ^ (-e^^(7(e~^^x,0,/i)) rfs 

T 

and this last limit is equal to 0, since the numerator is bounded in T for any given (x, n) G 
]^2p+i+s_ So we conclude that / G ker A, and thus that ker^ = range £^,1/- This establishes 
fl3.8|l , and concludes the proof of Proposition 13.71 I 
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